clear all;clc;

v00 = load('v_matrix.txt');
vx   = v00(:,1);
v0x  = v00(:,2);
vautx= v00(:,3);

q00 = load('q_paid_matrix.txt');
q_nodefx = q00(:,2);

na=21;
nb=61;

for i_b=1:nb;
for i_a=1:na;
        v(i_a,i_b)=vx((i_b-1)*na + i_a);
        v0(i_a,i_b)=v0x((i_b-1)*na + i_a);
        vaut(i_a,i_b)=vautx((i_b-1)*na + i_a);
        q_nodef(i_a,i_b)=q_nodefx((i_b-1)*na + i_a);
end;
end;

